Condensate Fraction of a Two-Dimensional Attractive Fermi Gas 
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We investigate the Bose-Einstein condensation of fermionic pairs in a two-dimensional uniform 
two-component Fermi superfluid obtaining an explicit formula for the condensate density as a func- 
tion of the chemical potential and the energy gap. By using the mean-field extended BCS theory, we 
analyze, as a function of the bound-state energy, the off-diagonal long-range order in the crossover 
from the Bardeen-Cooper-Schrieffer (BCS) state of weakly-bound Cooper pairs to the Bose-Einstein 
Condensate (BEC) of strongly-bound molecular dimers. 
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In the last few years several experimental groups 
have observed the crossover from the Bardeen-Cooper- 
Schrieffer (BCS) state of Cooper Fermi pairs to the Bose- 
Einstein condensate (BEC) of molecular dimers with 
ultra-cold two hyperfine component Fermi vapors of K 
atoms [H 0, Q and 6 Li atoms [3, A finite condensed 
fraction has been detected also in the BCS side of the 
stimulating a debate over its interpre- 
©]. Extended BCS (EBCS) equations 
have been used to reproduce density profiles 
13| and collective oscillations [lj] of these Fermi gases. 
Improvements based on Monte Carlo fitting and super- 
fluid dynamics have been considered [HI, [l6|, showing 
that EBCS is quite accurate. As this EBCS mean-field 
theory is defined for any value of the coupling, it provides 
an interpolation between the BCS weak-coupling regime 
and the BEC strong-coupling limit [l7|, EH ■ Despite well 
know limitations [TtJ the EBCS theory is considered a 
reliable approximation for studying the whole BCS-BEC 
crossover at zero temperature, giving a simple and co- 
herent description of the crossover in terms of fermionic 
variables. 

Recently, within the EBCS scheme, we have derived 
[l9| an explicit formula for the number of condensed 
fermionic pairs in the uniform BCS ground-state. We 
have used the EBCS equations to study the behavior 
of this condensate fraction as a function of the inter- 
atomic scattering length in the BCS-BEC crossover: from 
the BCS regime crossing the unitarity limit to the BEC 
regime. With no fitted parameters, we have found a re- 
markable agreement with experimental results Q indi- 
cating a relevant fraction of condensed pairs of 6 Li atoms 
also on the BCS side of the Feshbach resonance. Note 
that our analytic formula of the condensed fraction 
of the three-dimensional (3D) attractive Fermi gas has 
been obtained independently also by Ortiz and Dukel- 
sky [2(| and Ohashi and Griffin [21] . The formula has 
been then compared with Monte Carlo calculations by 
Astrakharchik et al. [22J: these Monte Carlo results 
show that the analytical formula slightly overestimates 
the condensed fraction of Fermi pairs. 

For repulsive Fermi gases, it has been predicted that 
a reduced dimensionality strongly modifies density pro- 



files [23j, [2J, 12a, 129], collective modes [27| and stability 
of mixtures [28i, 29] . In this paper we calculate the con- 
densate fraction of Fermi pairs in a strictly 2D attrac- 
tive Fermi system. It is well known that purely attrac- 
tive potentials have bound states in ID and 2D for any 
strength, contrary to the 3D case [3(3]. It follows that 
a bound state appears immediately as the two-body at- 
traction is introduced in a 2D Fermi gas [3ll [Hj]. As 
discussed by Marini, Pistolesi and Strinati [33|, for a 2D 
Fermi superfluid in the EBCS theory it is the value of the 
bound-state energy that determines the crossover from 
BCS state of weakly-bound Cooper pairs to the BEC of 
strongly-bound dimers. We shall show that the conden- 
sate fraction of Fermi pairs, which can be expressed in 
terms of elementary functions, increases smoothly during 
the 2D BCS-BEC crossover and becomes equal to 1/2 (all 
pairs are condensed) only for a very large bound-state en- 
ergy. 

The Hamiltonian density of a dilute and interacting 
two-spin-component Fermi gas in a box of volume V is 
given by 



H 



2m 



Y, ^V^+S^+^t . (1) 

'=t,l 



where Vv( r ) is the field operator that destroys a Fermion 
of spin a in the position r, while ip"^ (r) creates a Fermion 
of spin a in r. The attractive inter-atomic interaction 
is described by a contact pseudo-potential of strength g 
(g < 0). The number density operator is 



h(r) = Y t + ( r )i( r ) 

<r=T4 

and the average number of fermions reads 
N = d 3 r (n(r)) . 



(2) 



(3) 



The interacting term can be treated within the mean- 
field Hartree-Fock approximation, namely 



(4) 
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and the Hamiltonian density ([I]) is diagonalized by us- 
ing the following Bogoliubov representation of the field 
operator 
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(5) 



in terms of the anticommuting quasi-particle Bogoliubov 
operators bu_ a . The quasi-particle amplitudes Uk and v k , 
such that + v\ = 1, are obtained by imposing the 
minimization (30l | of the thermodynamic potential 



d 3 r (W(r) -/zn(r)) 



(6) 



and then subtracting this equation from the gap equation 
[ToL [Til [l2j , one obtains a regularized gap equation 



E 



2m 



(13) 



In the two-dimensional continuum limit ^ k — » 
y/(2vr) 2 /rf 2 k -» V/(27r) / fcdfc, taking into account the 
functional dependence (jlip of the amplitudes Ufe and Wfc 
on n and A, the Eq. (| 13|) gives 




(14) 



where fi is the chemical potential, fixed by the average 
number N of fermions. At zero-temperature the average 
of quasi-particle Bogoliubov operators is given by 



(7) 



where E k are the quasi-particle energies and 0(x) is 
the Heaviside step function. Neglecting the Fock terms 
{tp^tpa-i) with a ^ a' [13], after minimization of ((6|) one 
recovers 35, 36] the standard BCS equation for the num- 
ber of particles 
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and the familiar BCS gap equation 
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Here 



h 2 k 2 
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(8) 



(9) 



(10) 



while the number equation |8]) becomes 




(15) 



These two equations are exactly those found in the ap- 
pendix B of the paper of Marini, Pistolesi and Strinati 
[33| ]. We observe that, for a 2D inter-atomic potential 
described by a 2D circularly symmetric well of radius 
i?o and depth Uq, the bound-state energy es is given by 
e B ^ h 2 /(2mR%)ex-p(-2h 2 /(mU Rl)) with UqE 2 -► 

As previously stressed, several properties of ultra-cold 
Fermi gases have been investigated in the last few years 
by using the EBCS equations [I3,[l3|- Here we analyze 
the condensate fraction of fermionic pairs that is strictly 
related to the off-diagonal long-range order (ODLRO) 
[37l | of the system. As shown by Yang 38] , the BCS state 
guarantees the ODLRO of the Fermi gas, namely that, 
in the limit wherein both unprimed coordinates approach 
an infinite distance from the primed coordinates, the two- 
body density matrix factorizes as follows: 



and 



^+(^+(4)^(^(1-2)) 

= (^+(ri)^+(r 2 ))(^(r 1 )Vi T (r 2 )). 
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The chemical potential fi and the 



gap energy A are obtained by solving equations ([5]) and 
([5]). Unfortunately, in the continuum limit, due to the 
choice of a contact potential, the gap equation diverges 
in the ultraviolet. This divergence is logarithmic in two 
dimensions and linear in three dimensions. 

As discussed in [33[ , quite generally in two dimensions 
the bound-state energy €b exists for any value of the 
interaction strength g. For the contact potential the 
bound-state equation is 
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The largest eigenvalue A^o of the two-body density matrix 
(fTB"]) gives the number of Fermi pairs in the lowest state, 
i.e. the condensate number of Fermi pairs [TJ HH HI]- 
This number is given by 



(17) 



N = I d 3 n d 3 r 2 |(^( ri )V T (r 2 )>| 5 
and it is straightforward to show [39( that 



N = J2 u l y 



In the two-dimensional continuum limit we find 
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Finally, by using Eq. (fT5)) and Eq. (j2"0"]) we obtain a 
remarkably simple formula for the condensed fraction 
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This is the main result of the paper. Nicely, in Eq. (|2lj) 
the condensate fraction depends only on the parameter 
xq = /i/A. In the weakly-bound BCS regime (xq 3> 1) 
the condensed fraction no/n goes to zero, while in the 
strongly-bound BEC regime (xq <C —1) the condensed 
fraction no/n goes to 1/2, i.e. all the N/2 Fermi pairs 
belong to the Bose-Einstein condensate. 

In 2D the Fermi energy ep = h 2 k F /(2m) of a non- 
interacting Fermi gas is given by ep = irh 2 n/m. It is 
convenient to express all relevant energies in terms of 
the Fermi energy ep. In this way these scaled quantities 
depend only on the parameter Xq. In particular, we find 



a/1 + x 2 + x 



(22) 



A = 2 



+ , (23) 



and also 



cp 



= 2x 



U/l + *§-*„) . (24) 



All these quantities are parametrized by Xo- It is then 
quite easy to plot the scaled chemical potential fi/ep and 
the scaled energy gap A/ep as a function of the scaled 
bound-state energy es/^F- 

Fig. 1 shows that the scaled energy gap A/ep (solid 
line) grows by increasing £_b/£f while the scaled chemi- 
cal potential n/ep (dashed line) decreases. The chemical 
potential fi is zero when the bound-state energy es is 
equal to 2ep . For larger values of es the chemical po- 
tential fi becomes negative. Physically, we can say that 
the value 2ep discriminates between the BCS regime of 
weekly-bound fermionic pairs (0 < es < 2ep and u > 0) 
and the BEC regime of strongly-bound fermionic pairs 
(e B > 2e F and fi < 0). 

By using Eq. (j2T|) and Eq. |22|) we can also plot 
the condensate fraction no/n as a function of the scaled 
bound-state energy cb/^f- 

Fig. 2 shows the condensate fraction no/n of Fermi 
pairs. The fraction is zero when the bound-state en- 
ergy €b is zero. For small values of eb/ef the con- 
densed fraction has a very fast grow but then it reaches 
the asymptotic value 1/2 very slowly. Note that at 
£s/cf = 2 (where fi = Xo = 0) the condensate fraction 
is tt/8 ~ 0.39. 
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FIG. 1: Energy gap A (solid line) and chemical potential fi 
(dashed line) in the uniform two-component dilute 2D Fermi 
2clS ELS Qb function of scaled bound-state energy es/eF. The 
horizontal dotted line simply shows the zero. 




FIG. 2: Condensate fraction no/n of Fermi pairs (solid line) in 
the uniform two-component dilute 2D Fermi g function 
of scaled bound-state energy es/^F- The horizontal dotted 
line shows the asymptotic value no/n = 1/2. 



To the sake of completeness, we observe that in the 3D 
case the condensate fraction can be expressed in terms 
of the dimensionless inverse interaction parameter y — 
(^af)" 1 , where kp is the Fermi wave vector of 3D non- 
interacting fermions and of is the 3D s-wave scattering 
length of inter-atomic potential 19]. In this case, y G 
(—oo, +oo) and only for y > there is the formation of 
a dimer with bound-state energy eb — h 2 /(2map). The 
condensed fraction goes to zero for y — > — oo and to one- 
half for y +oo [191 ]. 

In conclusion, by using the mean-field extended BCS 
theory and the concept of off-diagonal long-range order, 
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that is the existence of a macroscopic eigenvalue of the 
two-body density matrix, we have obtained a remarkably 
simple formula for the condensate fraction of fermionic 
pairs in a uniform 2D Fermi gas. Contrary to the 3D case, 
in the 2D case a bound state appears immediately as the 
two-body attraction is introduced. As a consequence, the 
crossover from the BCS state of weakly-bound Cooper 
pairs to the BEC of strongly-bound dimers is induced 
by the increasing of the bound-state energy. We have 
show that the condensate fraction of Fermi pairs grows 
smoothly during the 2D BCS-BEC crossover, but only for 
a very large bound-state energy one gets a quasi complete 
condensation. It is important to stress that our predic- 
tions on the behavior of the condensed fraction in a 2D 



attractive Fermi gas can be surely compared with Monte 
Carlo calculations, as done in the 3D case (22J. On the 
other hand, it could be difficult to compare the theory 
with experiments. In fact, strictly 2D superfluid Fermi 
gases have not yet been achieved: a 2D configuration re- 
quires na 2 H <C 1, where n is the 2D number density and 
an is the characteristic length of a very strong harmonic 
confinement along one of the three axes [4(| ■ In addition, 
in two-dimensions it could be more problematic than in 
three-dimensions to modify and trigger the bound-state 
energy [4l| . Surely in the next future these issues will be 
faced and probably the obstacles will be overcomed. 

The author thanks A. Parola and F. Toigo for useful 
discussions. 
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